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The present paper studies the symmeries of the Hubbard model of electrons with generally n-fold 
orbital degeneracy. It's shown SUd{2n) and SUc{2n) symmetries hold respectively for the model with 
completely repulsive or attractive on-site interaction and that with partly attractive interactions. 
An extended Lieb-Mattis transformation is given to map these two symmetries into each other. 
The sub-symmetry 5(/^">(n) (g) SU^/\n) is found to be shared by the two models with arbitrary 
chemical potential ^. By assuming at most two electrons on each site it's found that SUd{2n) and 
SUc{2n) both exist in each kind of the two models and consequently lead to a larger symmetry 
SUd{2n) X SUc(2n). Another underlying symmetry (SC/i^' (2)-p x ... x SU^''\2)r) 8 {SU!,''\2)v x 
... x SUi°\2)v^ is also revealed for the unified U model under the excluding. The symmetry is valid 
for the partially attractive model with chemical potential fi — —U. 
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I. INTRODUCTION 



Recently considerable attention has been directed to the studies on correlated electrons in the present of orbital 
degree of freedom which is relevant to transitional-metal oxides Cgo materials Q and artificial quantum dot 

arrays 0. Apart from the numerical ||] and perturbative |^ works theories based on symmetries were presented 
for one-dimensional models of these systems. An SU{4) theory describing spin systems with orbital degeneracy 
was proposed (l^Jll[l for a theoretical understanding of the observed unusual properties. The ground-state phase 
diagrams for the system with a symmetry breaking of SU (4) — > SU{2) x SU{2) were discussed in ref. [Hl|. For the 
2-fold orbital degenerate Hubbard model a recent paper [n3| presented the SU (4) theory and showed the underlying 
SUd{^) symmetry of spin-orbital double and a charge SUc{4:) symmetry with an extended Lieb-Mattis transformation 
mapping those two 5^7(4) generators into each other. On the basis of elementary degenerate perturbative theory, it 
was also shown that the effective Hamiltonian is equivalent to the 5*0(6) and SU{4:) Heisenberg models respectively 
at half- filling and quarter- filling with strong coupling. In ref. jl4j, the one-dimensional S'[/(4) Hubbard model is 
extensively studied on the basis of Bethe ansatz solution. As for the symmetry theory of one-dimensional Hubbard 
model without orbital degeneracy, it has been well-investigated. Yang introduced the pairing operators and so 
constructed the symmetry SU{2) x SU(2) Based on the symmetry Fabian H.L.E/31er et. al. discussed 

the completeness of the Bethe ansatz solutions, M.Pernici showed the off-diagonal long-range order, D.B.Uglov 
and V.E.Korepin presented the Yangian symmetry Y{sl{2)) ®Y{sl{2)) |18|, and Fabian H.L.E/31er and Holger Frahm 
considered the density correlations [19| . And it has been argued that the two-dimensional single-band Hubbard model 
has approximate 5*0(5) symmetry pffl. But the research works on the Hubbard model with orbital degeneracy are 
still in accumulation. In the present paper we study the symmetries of the Hubbard model of electrons with generally 
n-fold orbital degeneracy. We show and clarify that the SUd{2n) and SUc{2n) symmetries hold respectively, unlike in 
the simple Hubbard model for which both of the two symmetries are valid, for model with unified on-site interaction 
and that with partly attractive interactions. But the sub-symmetry SUjf^ (n) (g) SUjf^ (n) is found both valid for the 
two models and for arbitrary chemical potential /i. An extended Lieb-Mattis transformation as in is given to 
map these two symmetries into each other. By assuming at most two electrons on each site, on the basis of which 
the Bethe ansatz can be applied in one-dimensional model, we find the SUd{2n) and SUc{2n) symmetries both exist 
in each kind of the two models so we have a larger symmetry SUd{2n) x SUc{2n). Under the exclusion another 

underlying symmetry (SU^''\2)-p x ... x SU^''\2)-p) (g) (5'J7i°''(2)p x ... x SU^°\2)-p), which is not included in the 
SUd{2n) X SUc{2n), is also revealed for the unified U model with = U so that the model totally possesses the 
symmetry SUd{2n) x SUc{2n) x [(S't/i''' (2)p x ... x SU^''\2)v) ® (S'C/i°\2)p x ... x SU^°\2)-p)]. The underlying 
symmetry is also valid for the partially attractive model with ^ = —U. 
Consider the n-fold orbital degenerate electrons with states 
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|1)=|1,T), |2)-|i,i), 

. . . , 

|2n-l) = |7^,t), |2n) = |n,i), (1) 

where in a state \l,<y) I denotes the Zth orbital component and a | label the spin components. We start with a 
general Hamiltonian with n-fold orbital degeneracy expressed by 

x.x' a X a^a' x^a 

where C^{x) creates a fermion of state \a) at site x and na{x) is the corresponding particle number operator. The 
notation of site is not restricted to one-dimensional case. 

II. SUd{2N) and SUc{2N) SYMMETRIES 

Besides the U (1) symmetry there exist two kinds of SU(2n) symmetries for the orbital degenerate Hubbard model. 
We define 

X 

Dm = -Nra+i, iV™ = ^ C+ (x) Cm (x) , (3) 

X 

they fulfill the commutation relations 

[Ess',Ett'] — Ss'.tEst' — Ss,t'Ets', 

[Dm, Ess'] — {Sm,s ~ Sm.s' ~ <^m+l,s + Sm+l,s')Ess' ■ (4) 

These operators can construct an SU{2n) Lie algebra 

SUd{2n) : {Dm, Ess' \ m = 1, ...,2n - 1;1 < s ^ s' < 2n} , (5) 

with n{2n — 1) of Egg' and there are totally (2n)^ — 1 of generators. {Dm} forms the commuting Cardan subalgebra 
of rank 2n — 1. The i?g(s_|_i)'s are the generators related to the simple roots. For any of txx' and fj,, all generators of 
SUd(2n) commute with the Hamiltonian (|^) with imified on-site interaction Uaa' — U , so we have SUd{2n) symmetry 
of spin-orbital double in this case. 
Let us define another set of operators 

X 

Pa2k ^^f{x)C2k{x)C2k+l{x), 

X 

0" = \ Y}^Ux)Cm{x) + C+^,{x)Cm+l{x) - 1], (6) 

X 

where fc = 1, n, m = 1, 2n — 1, f{x)'^ = 1, and f(x + S) = —f{x) for any site x and nearest-neighbor x + 6. The 
above operators can realize another SU{2n) Lie algebra which we shall denote by SUc{2n): 

{Qm,F^,F^^}. (7) 

{Qm} is the Cardan subalgebra. The -Fa^^.^jand Fa^^ are the generators related to the simple roots, other generators 
relating to positive roots can be obtained by Fa^+aj = [Fai,Faj], the generators with negative roots will be F^a = 
{Fa)'' . If we assume that the on-site coupling Uaa' — U for the states labeled by a, a' with different spin components 
while Uaa' — —U for states with the same spin components 

Uaa' = U for odd-even pair a, a', 

Uaa' = ^U for odd-odd or even-even pair a, a', (8) 
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and the amplitudes t^x' of odd-neighbor hopping are real and those of even-neighbor hopping are imaginary 

t*^, — txx' when x-x' is odd neighbor, 

= —ixx' when x-x' is even neighbor, (9) 

we will have the following relations 

[i7,^^„,J = 2(M-C/)F„,,, [H,F^,,_,]^-2{^i-U)F^,,_,. (10) 

If the chemical potential n ~ U, the Hamiltonian will commute with all the generators of SUc{2n) so the model has 
the charge SUc{2n) symmetry. In terms of the partially attractive (|^), it can be easily proved that such a Hamiltonian 
with fi = U has the half-filled form 

H =-J2Y.{^xx'C:{x)Ca{x')+t:^,C+{x')Ca{x)) Uaa'{na{x) - - 1). (11) 

x^x' a X a^a' 

The usual Hubbard model with the nearest-neighbor hopping is included in the class of (||). Whether f{x) can be 
well defined depends on the lattice structures, if in a lattice any sum of the nearest-neighbor pair 6 + 6' is not a 
nearest- neighbor and there are even sites on every perpendicular directions, f{x) can be well defined as 

/(x) = exp(i7r • x), tt = (tt, tt, ...), x = (a;i, a;2, ...), (12) 

Xi is the ith components of the site coordinate in the lattice basis. Such lattices (bipartite lattice) can be simple 
squared, centered squared in 2-dimension; simple cubic and body-centered cubic in 3-dimension. For 1-d case and 
even total sites, f{x) = exp(j7rx). 

It should be noted that the SUd(2n) does not commute with the partially attractive Hamiltonian which has only 
the SUci^n) symmetry, neither does SUd^n) with the unified U model which possesses the SUd{2n) symmetry. The 
two kinds of symmetries can be mapped into each other by an extended Lieb-Mattis transformation 

Ca{x) 1-^ exp(i7r • x)Cg^(a;) for even a, 

Ca{x) 1-^ Ca{x) for odd a. (13) 

The transformation leaves the hopping term (^ invariant and changes the sign of Uaa' with odd-even pair a,a' in ( [ll| ) 
so maps the model ( |ll[ ) into the unified Uaa' — —U model. 

The particle number of each states and the total spin can be expressed by 

2n-l 
m — 1 

2n-l 2n-l 

^1 = + - H mDm)l2n, j < 2n, (14) 

ni—j m—1 



Stotal = ^2m-l/2, (15) 

m—1 

where Ng is the number of total electrons. 

III. su'-^\n) (g) su';^\n) symmetry 

Unlike in the transitional Hubbard model for which both of the two symmetries are valid for on-site attractive 
and repulsive interactions, as we can see in previous section, for the orbital-degenerate Hubbard model the SUd(2n) 
symmetry of spin-orbital double only holds for the unified U Hubbard model with arbitrary chemical potential ji 
whereas the charge SUc{2n) symmetry merely exists in the other partly-attractive half-filled model of which the 
chemical potential is /i — U. But considering that SUd(2n) and SUc{2n) share some common generators E2u^2u> and 

E2i^-i,2u' we will find the shared sub-symmetry Sul[\n) (g) Sujf\n) yields for the two models with 
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su. 



SU^/^ (n) : {I?2.,2.+2, £;2.',2." W<n-l,l<:y' <n}, 



(16) 



where Djn.m+2 — Nm — Nm+2- Especially, we shall illustrate that the symmetry is valid for the partially attractive 
model with arbitrary chemical potential /i. 

As the hopping term in the Hamiltonian and the chemical potential term M X^a; a SUd{'2n) invariant, 
whether the Hubbard model possesses SUd(2n) symmetry depends on the commutation relation of the on-site inter- 
acting term and the SUd{'2n) generators 



E^,^ Uaa'naix)na'ix) ^ C:^{x)C„i+k{x)(^ ^ Hg' {x)Um+k,a' - na>{x)Um,a') 

a'^m--\-k a' 

+ ( ^ na{x)Ua,m+k - na{x)Ua^rn^C:^,ix)C.m+kix) 



(17) 



ay^7n-\-k 



where 



Et^Y.^+{^)^m+k{x). 



Surely for the case of unified Uaa' — U we easily find that the above commutation vanishes so that we have the 
SUd{'2,n) symmetry in this case, as is obtained in the second section. Although for the partially attractive model the 
above commutator does not go null for all fc's and consequently we do not have SUd{2n) symmetry, the case with 
even k will be an exception. From the partially attractive (^) we find for even fc's 

Ua^m+k — Ua,ni — 



,m-\-k ''^771+ k{x)Ura+k ,1 

a^m-\-k a^m 

= {nm{x) ~ nm+k{x)){~U). 

Then eq.(|l7|) becomes 

Emi X! Uaa'na{x)na'{x) 



= C+(x)C„+fc(a;)(n„(a;) - n„+fc(a;)) (-J7) + (n™(x) - nyn+k{x))C^{x)Cra+k{x){~U) 
= {-C+{x)C„,+u{x) + C+{x)C„,+k{x)){-U) = 



where we have used relations: rimCr, 



0, TijYiC^ — Crn'^m 



Cm- with even fc's and odd (or 



even) m's correspond to the SU^"'{n) (or SU^^'''{n)) generators. As a result, the SU^f\n) eg) SU^f (n) generators in 
( |l6| ) commute with the Hamiltonian of the partially attractive model for any chemical potential fi. Therefore the 



SUjf\n) (X> SU^f{n) symmetry is shared by the two models both for arbitrary fi. 



t(o) 



IV. AT MOST TWO ELECTRONS ON EACH SITE 



The application of Bethe ansatz method to the 1-dimensional degenerate Hubbard model is based on such an 
assumption that prevents scattering process involving three or more electrons on one site 
Hubbard model the configurations of more than two electrons on one site are excluded automatically by the Pauli 
principle. In the continuum limit and for small densities or [/ 3> i in the lattice model, the unwanted configurations 
in degenerate Hubbard model become negligible, so the Hamiltonian with three-electron configurations excluded will 
describe the system well. If we exclude more than two electrons on each site, we will find that the Hamiltonian has 
both SUdi^n) and SUc{2n) symmetries and furthermore a larger symmetry SUd{'2n) x SUc(2n). In addition, we will 

find an underlying symmetry (S-f/i"^ (2)p x ... x SU^''\2)v) ® {SU^°\2)v x ... x SU^°\2)v). 
Consider Uaa' — U case, the Hamiltonian reads 




. For the traditional 
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n=-tY. ^PC+(x)a(x')P + i7^ ^ n„(x)n„,(x) (18) 

<x,x'> a X a^a' x,a 

where < x, > represents the nearest neighbor sites and let us define 

•^c,,., =J2exp{tn ■ ^)VC+,_,{x)C+,{x)V, (19) 

X 

•^a2. = J2 ^^P(*'^ • ^)^<^2fe {x)C2k+l {x)r, (20) 

X 

where the operator V projects onto the subspace of states having at most two electrons on each site Other 
generators can be given from the ^aat-i •^a2k the SUd'in) in section II. The V operator excludes such 

terms as nsC^,C^,,, CgCs'Tis" with three different s,s' and s" so that 

[VCt{x)C+ix)r,J2Y. Mx')na'{x')] = ^2VCt{x)Ct,{x)V, (21) 

x' a^a' 

thus we have the relations similar to (|l^) 

[n,To.,,]^2{^l-U)T„,,, [n,T^,,_,] = -2(n-U)To,,,_,. (22) 

Set n = U and the ^a^k^i ^^'^ ^a^k ^^^^ commute with the Ji. J-a^k-n ^oi2k ^'^d Ns — Ns' can generate the SUc{2n) 
symmetry so that TC has both SUd{2n) and SUc{2n) symmetries, 

SUc{2n) : {Qrn, , (23) 

SUdi2n) : {Vm, £ss'} , (24) 

where Qm = VQmV and the SUd{2n) generators are correspondingly revised to be 

Ess' = VEss'V, = VD^V. (25) 

It also can be similarly shown that both of the two symmetries hold for the "P-modified model with partially attractive 
Uaa' — U, — U. As the two symmetries both hold for each of the models we can construct the larger symmetry 
SUd{2n) X SUc{2n) for each of them. If we do not exclude more than two electrons on each site, these two symmetries 
respectively belong to different models. 

Besides SUd{2n) x SUc{2n) there exist some less obvious symmetries. Define 

HtL' = E exp(»'r ■ ^)PC+,ix)C+,, ix)V, = (-^^fc.LO^ (26) 



•^2fe-l,2fe'-l = EexP(*'^ ■ x)7'C+_i(x)C2+,_i(x)7', = (^i^!i,2fc._i)^ (27) 

X 

such generators are not included in ( p3| ) which contains pairing operators ^^^^only with odd-even pair ss' . It can be 
easily verified that these operators also commute with the unified \J model ( |l8|) with [i~U so we find an underlying 

symmetry {SVf^{2)v x ... x Svf^{2)v) ® {SVf{2)v x ... x Su'f\2)v) with 

SUf\2)v:[Q2K'iU', ^Sfe" -^Sfc'} fc^fc'' Kk'^\.....n, (28) 



5C/(°)(2)p:{Q2fe_i,2fc,_i, ^i+\2/c'-i, -^iM'-i} Mfc', fc,fc' = l,...,n. (29) 

There are respectively C^^ ~ n(n — l)/2 of the SUc'^\2)-p and SUc°\2)-p symmetries. The extended Lieb-Mattis 
transformation (|l^) maps the above symmetry into itself. An revised Lieb-Mattis transformation mapping into 
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the corresponding SUji'\2)-p and Sujf\2)-p will involve a third kind of Hamiltonians with partially attractive 
Uaa', which differs from what we discussed before. So finally we have the symmetry SUd{2n) x SUc{2n) x 

{SU^''\2)'p X ... X SU^''\2)'p) ® {SU^°\2)-p x ... x SU^°\2)v)\ for (|l|) with fi^U. But for Uaa' = U, ~U case the 

underlying symmetry is valid for another chemical potential. The Hamiltonian under that exclusion is 

<x,x'> a X a^a' x,a 

where the partially attractive interaction Uaa' is also defined by (||) . Compared with ([lo| ) the commutation relations 
are different 

[W':^£L]--2(M + ^)^i:,L> (31) 



[^',^L2fe']=2(A* + C/).Fy2,„. (32) 



The different sign of U comes from U2k,2k' = —U while in eq.(fO) it's U2k,2k±i = U. Therefore the symmetry 
{SU^''\2)v X ... x SU^''\2)v) ® (5[/i°^(2)p x ... x SU^°\2)v) holds for ^ = -U. It should be noted the symmetry is 
valid under the exclusion, without the T-'-exclusion its generators will commute with neither of the two kinds of models. 
And unhke SU^°\n) (g, SU^''\n) C SUdi2n) in section III, none of the SU^°\2)v or 5f7^^(2) is any sub-symmetry of 
SUc{2n). 

Considering the commutations J-a2k ^^'^ ^ss' , we have 

[SU^/''\2),SU^/''\2)] = (33) 



[SU^/'''\2),SU^'"'"'H2)] = 0, for seperate pairs (ss') and (s"s"'), (34) 

but 

[SU^/''\2),SUj:'''"\2)] / 0, 

where SUjf \2) and S't/i'** ^(2) are sub-symmetries involving only the states s and s' . The whole symmetry cannot 
be written in a direct product SUd{2n) ® SUc{2n) but SUd{2n) x SUc{2n). This is also a difference from the 
single band Hubbard model of which the symmetry in our notation is a direct product of the two S'[/(2)'s: i.e., 
50(4) ~ SUd{2) ® SUc{2). Therefore for the one-dimensional model ( p^ ) which can be solved by the Bethe ansatz 
p2| , the discussion on the completeness of the Bethe ansatz solution and the off-diagonal long-range order will be 
expected quite different since in the single band case it's based on the vanishing commuation of the SUd{2) and 
5(7,(2). 



V. BRIEF SUMMARY 



In summary, we studied the symmetries of the Hubbard model of n-fold orbital degenerate electrons. We show and 
clarify that the SUd{2n) and SUc(2n) symmetries hold respectively for the model with unified on-site interaction and 
that with partly attractive interactions. An extended Lieb-Mattis transformation is given to map these two symmetries 
into each other. But the sub-symmetry SUjf\n) (g) SUj^°\n) is found to be possessed by the two models and both 
for arbitrary chemical potential fi. By excluding more than two electrons on the same sites we find the SUd{2n) and 
SUc{2n) symmetries both exist in each kind of the two models, so we have an enlarged symmetry SUd{2n) x SUc{2n). 
Under this exclusion, another underlying symmetry {SU'f\2)-p x ... x SU'f\2)-p) ® {SU'f\2)-p x ... x SUc°\2)j>) is 
also found for the unified U model with chemical potential fi — U, and consequently this model has the symmetry 

SUdi2n) X 5;7c(2n) x ["(5C/i''^ (2)p x ... x SU^''\2)-p) (g) iSU^°\2)-p x ... x SU^°\2)v)\. The underlying symmetry is 
valid for the partially attractive model with chemical potential jj, — —U. 
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